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The relaxation of the harmonical oscillator, being in contact with the ther- 
mostat whose excitations can be treated as bosons is studied. Temperatures 
CN , of the oscillator and thermostat are supposed to be different at t=0. An 

exact temporary behavior of the inverse oscillator's temperature is obtained. 
In low interaction limit the temperature behavior is exponential. When the 
interaction is increasing the temporary behavior display oscillations. Narrow 
peak of the oscillator susceptibility is discovered. The exact solution obtained 

'& . is compared with Bloch-Redfild approximation, 
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If the reader knows some publications related to the susceptibility peak described below, 
please, send me the corresponding reference. 

The problems of kinetics are known to be one of the most complicated problems of 
physics. Existing approximate methods usually do not allow one to estimate the quality 
of the approximation, and the plausibility of the results is the main criteria of the validity 
^ ■ of the approximation made. For this reason the simple exact-solvable models of dynamical 
system and thermostat are of particular interest. The exact solution can be compared in 
detailes with the solution, obtained by means of the approximate method. This allow one to 
evaluate used approximation in a correct way. Moreover the exact solution can demonstrate 
some new unexpected phenomena. We consider the resonance of the susceptibility described 
below to be one of them. In this paper we describe the model of this kind with the harmonical 
oscillator as a dynamical system and system having boson -like excitations as a thermostat. 
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I. THE PROBLEM DESCRIPTION AND RESULTS 



This model was considered by Bogolubov |1| for classical vibrator. The typical example of 
such thermostat is the phonon system of the crystal. The model's Hamiltonian in frequency 
units has form: 

H = H p + H s + H sp (1) 

ft ' 

N 

H P = Y, u q A + q A q , H s = lo s B+B 

9=1 

N b 

9=1 
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Here H p - the Hamiltonian of the bosonic thermostat, u q -the boson's dispersion (for 
the phonon thermostat - the phonon frequencies) H s - The Hamiltonian of the harmonical 
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oscillator, H sp — the Hamiltonian of the interaction , N- number of bosonic modes, Cq - 
the interaction constant. The inverse proportionality of this constant to the number of 
bosonic modes N is typical for all known interactions. A't, A q (B + , B) -the operators of 
creation and annihilation of bosones in q-mode (oscillator's energy quantums). The following 
commutation rules are taking place: 

[A q ,A+] = 5 q , q ,, [B,B+] = 1 (2) 

We accept the simplest form for the bosonic spectrum: 

2W 

u q = u s + —(q-N/2-3/4), q = l,...,N (3) 

So, the oscillator frequency is coincide with bosones spectrum center. Suppose the oscillator 
and thermostat has at t — the temperatures T p q and T s q respectively. So, the problem is 
to solve the Lioville equation for the density matrix p: 

ip = [H, p] (4) 

with the initial condition: 

p(0) = ^ 1 exp(-^±^-^), (5) 

where 

Z = Spexp(-5i±^-|l) (6) 

Here we consider the temperature of interaction equal to that of thermostat. The unequi- 
librium state of this type can be obtained in a following way. Let the oscillator's frequency 
depends linearly on some external parameter h: u s = uh. Let this parameter has the 
value hi for a long time and our system reached the equilibrium, having the density matrix: 
~ exp — [u' s B + B + H sp + Hp\/TpQ, oo' s = Qh! . Now let this parameter change instantaneously: 
h! — > h, and h — h! = 5h and, therefore, u' s — » u s . After that the Hamiltonian of the system 
has the form (1). It is easy to see that right after jump of the parameter h, the density matrix 
can be represented in the form (5), with T s0 = (u s /u' s )T p0 . Thermostat is supposed to be 
microscopic, so all calculations below we perform in thermodynamic limit N — > oo. Lastly 
note that the above model can describe the spin relaxation at temperatures low enough as 
compared with spin levels splitting in magnetic field. In this case the parameter h can be 
associated with magnetic field. 

Thus problem is described, we now itemize the main points of the solution. For this 
reason we make the following general remark. Consider the set of operators Oi , % — 1, ...M 
to be closed with respect to the commutation with Hamiltonian in the following sense: 

M 

[H,O t ] = J2^ k O k , (7) 
k=i 
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where - numbers. Take now the density matrix in the form: 

M 

a(t) = J2mO l , (8) 
i=i 

substitute it in (4). Making use of (7), we calculate commutators and equalize the factors 
corresponding to the similar operators. We obtain the following set of equations for Pi(t) 
functions: 

M 

z& = z' = l,..,M (9) 

fc=i 

After solving this set of equations we obtain the solution of Lioville equation in the form 
(8). Note now that the set of operators A q A q , A q B + , A q B, q = 1, N and BB + is closed 
in the above sense with respect to the commutation with Hamiltonian (1) and we can use 
the above method as follows. Let us put into consideration the auxiliary density matrix a: 




Substituting it in (4) and equalizing the factors corresponding to the similar operators 
we obtain the following set of equations: 

4 = C (7 9 - 7,*) (11) 

N 

x/j^-CbXfo-Tj) 

q=l 

liq = (UJ S - U q )j q + C (P q ~ p s ) 
l % = ~(U S ~ Uqhg ~ C (f3 q ~ (3 S ) 

We solve this system with initial conditions: 

(3 q (0) = u q /T p0 , (3 s (0)=uj s /T s0 , lq (0) = 1 ;(0) = C /T p0 (12) 

After that we use known property of Lioville equation (p. 170) : if matrix a satisfy the 
Lioville's equation, then the arbitrary analitical function of this matrix satisfy the Lioville's 
equation too. Taking this into account, it is easy to see that matrix p = Z^ 1 exp(— a) 
satisfy the Lioville's equation with desired initial conditions (5). The final solution of the 
problem we obtain by turning to the thermodynamic limit N —>■ oo. It is seen that functions 
P q , P s can be called the " dimensionless inverse temperatures" of the q bosonic mode and 
dynamical system respectively. We call ^ q functions the "dimensionless inverse temperatures 
of the interaction" with the bosonic mode q. An exact solution of (11) performed in the 
second section lead to the following formula for the temporary dependance of the inverse 
temperature of the oscillator T^ 1 = p s (t)/u s : 
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1 — cos Vlt 



n i + (nr-A(n)J 



(13) 



^ n l2 



(14) 



(15) 



Here p = 1/2VF - is normalised to unit density of bosons modes. If b/W << 1 one can 
extend the integration up to the infinity and drop the logarithmic term A(Q) in denominator. 
After that the integral can be calculated and we obtain: 
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— + — - — 1 - exp 



L p0 



1 



(16) 



The temporary dependances, calculated by formula (13) for different b/W , are presented 
on fig.l. It is seen that when interaction is small enough, the kinetics is exponential. In- 
creasing of the interaction is accompanied by occurrence of oscillatory kinetics. We present 
here the expression for the inverse temperature of the bosonic q- mode q 1/T g = I3 q /u) q : 



1 

a 



T, 



(17) 



1 



9 r W 
7T JO 
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T p0 T s0 J\ Nu q JUJo (fi 2 -cc; 2 )[l + (fir- A(fi)) 2 ] +(1 COS ^ 1 + (our - A(cu)) 2 

u 2 



(u q - uj s f + 2C 2 
If b/W « 1 this expression can be simplified: 
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[Nu q {l+u 2 T 2 ) 



r(l — cos ujt) + 
Formula (14) can be rewritten in the form: 

T - x = 2irbp (uj s ) 



sin uj \t\ 
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^1 — exp — 


T 


)] 



(19) 



Po(u s ) - density of the resonance bosonic modes. 

Let us now return to the dropped logarithmic term which can be of particular interest 
for the following reason. The above problem describes the 'system's response' on the jump 
of the parameter h . If we consider the 'system's response' as S(l/T s ) = (1/T S — 1/T p0 ), 
then S(l/T s ) = $(t)5h, where $(£) can be derived from (13). Let parameter h gains small 
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harmonical increment 5h LU exptujt. The corresponding small modulation S(l/T s ) w exp iujt 
can be expressed in terms of susceptibility: 

On the other hand, the linear response theory []2| (. 88) lead to the following relationship 
between the susceptibility and 'system's response' on the jump $(£): 

2ir J-oo nj + Jo 

The imaginary part of susceptibility can be explicitly calculated. It differs from zero only 
when —W<uj<W and is equal to: 

UJT CJ S 

Im Xu = Xoo ——, * , \ \ 9 Xc 



1 + (ujt - A(uj)) 2 T p0 cj s 

It is seen that susceptibility Im^^ exhibit two extraordinary extremums whose positions 
±a;o (when b/W << 1) can be defined from the equation: 

UJ T = A(Uq) 

It is easy to see that ujo ~ W — 0. So, if uj = ujq the modulation of oscillator's inverse 
temperature in the orthogonal phase (with respect to 8h), is Wt = (W/b) 2 /n » 1 times 
greater than the adiabatic modulation Xoo$h . If interaction is small enough, it is possible to 
estimate the absorption under the action of modulation Sh = Sh u sin out. The rate of energy 
changing dE/dt for whole system (oscillator+thermostat) can be calculated as: 

= ^Sp pH = Sp P~jj.Hs = u s Sh Sp pB + B = uJ s 5h < n > 

If the interaction is small, terms ~ 7 g can be omitted in p ~ exp(— a). In this case 
the average number of quantums < n > is defined by the Bose-Einstein function < n > = 
l/(e^ s — 1) with p s = uj s /T s depending on time. Linear response 8f3 s related to the modulation 
5h can be determined by above expression for susceptibility. Furthermore dE/dt should be 
averaged over the period of oscillations. For this reason only the component of 5(1/T S ) having 
the same phase as 5 h is important. It is defined by above expression for the imaginary part 
of the susceptibility. Taking all this into account it is easy to see that l^dE/dt^j is differes 
from zero only if < uj < W and in this frequency interval can be calculated as: 

/dE\_5ul^P_ U 2 T R = ^L x =-x h 

\ dt I 2T p0 (e? - l) 2 1 + [lot - A(u,)] 2 P " V 3 ~ Us " 

The frequency behaviour of the absorption is presented diagrammatically on fig. 1 (c). 
It is seen that the absorption exhibit narrow maximum at uj = ujq. The occurrence of this 
maximum is somewhat unexpected. This peak do not depend on details of the bosonic 
spectrum. 
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II. THE EQUATION FOR THE INVERSE TEMPERATURES 



In this section we solve the set of equations (11) which can be rewritten as: 

i q = _[(^ _ Ug f + 2C 2 ]I q - 2C 2 Q £ J„ q = l,...,N (20) 

ql=l 

L q = Iq ~~ lq 

As a starting point we find the solutions of this equation having the harmonical temporary 
dependance - normal modes exp(tQt), q — 1, N. 

After substituting in (20) we obtain the following equations for the amplitudes : 

e n = 2C ° S n = l N (21) 

where 

N 



S = E^ (22) 



q=l 

Substituting (21) in (22) we obtain for squares of eigen frequences Q 2 : 



N 9^2 

For the bosonic spectrum in the form (3) this equation have N undegenerated roots. 
Using (21) we obtain for the normal modes the following expression: 



e n = MS m 

q {uj s - u q ) 2 + 2Cq - VL 2 

Where N (Q) is defined as: 



1 



(uj s - u q y + 2c 2 - w 



2 



(25) 



(24) is the eigen mode of (20) only if f2 2 is coincide with one of the roots of (23). The 
initial conditions for quantities I q can be obtained from (12): 

7,(0) = 0, / ff (0) = 2iC7oa;.(^--^-) (26) 

Introducing the solution of (20) in the form of expansion over the normal modes, we can 
write down the solution of (20) as: 

/ 1 1 \ ^ e^. sin Qt 
I q (t) = 2zC u s (—- — £ (27) 



T p o T s0 J Qqf il 
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il runs over the set of eigen frequences obtained from (23). From (23) and (24) one can 
obtain: 



' 2C 2 



(28) 



Using this relationship we obtain the solution of (20) in the form: 

u„ ( 1 1 \ ^ N 2 (tt) sin tit 



L = -i- 



E 



C \T p0 T s0 J ^Q[{uj q -uj s ) 2 + 2C 2 ~Q 2 } 



(29) 



Taking into account (28), (27), (11) and initial conditions for l/T s = f3 s /uj s it is easy to 
obtain: 
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T s qJ 2C VTpo T s0 J n 



E^) 



1 — cos Qt 

n 2 



(30) 



The sum in this formula can be replaced by integral if density of squares of eigen fre- 
quences is known. We calculate this density as follows. The above eigen frequences arises 
in the following problem for eigen values: Me = f2 2 e, where matrix M is defined as: 



d r . 



V 



M = M + V 

o \ 



(31) 



d N J 



V = 2C 2 



W 2 



1\ 



(uj s - u q f + 2Cl - — (m + 1/2) 2 + 2C 2 



The diagonal elements d m are written for the bosonic spectrum in the form (3). Let us 
put into consideration the Green's function: 



G{E) = (E-M)- 1 



(32) 



Then the density of squares of eigen frequences of interest D(fl 2 ) is defined by known 
formula: 



D(E) = Im SpG(E), 



E = tt 2 + t5 



+0 



(33) 



Considering V as perturbation for M , one can write down the Dyson's expansion for 
Green's function: 

G = G + G VG + G VG VG + ... G = (E - M ) _1 (34) 

This series can be exactly summed and for the Sp G of interest we obtain: 

T(E) dT/dE 



SpG(E) 



2Cl \-T{E) 
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(35) 



where 



T(E) = 2Cl tt^t = 2C 2 Sp Go (36) 



For the bosonic spectrum in the form (3) in the thermodynamic limit one can write down 
the explicit expression for T(E): 

(2b 2 \ N 1 
r(E) = { "iV J JJ E - (W/N)\m + 1/2) 2 - 2b 2 /N = (3?) 

2 b *r ^=*^ +w 



lo E-(W£) 2 W\[E VE-W 
Factor (25) can be presented in terms of T(E) as: 



^ (n) = -A E = a2 (38) 



And taking into account (33), (35) we obtain the following relationship: 

dE 1 - cos \J~Et 



(39) 



Making use of (30), (39) it is now easy to obtain (13). Formula (17) can be obtained in 
the same manner. 



III. THE COMPARISON WITH BLOCH-REDFILD'S THEORY. 



This theory |2j (. 60) lead to the following kinetic equation for the diagonal elements of 
density matrix of above oscillator: 



Pn — Pn+lW / n+l-+n + Pn~\ — 1— »n Pn 



where W, 



n— mil • 



1 



n + 1 



•nil 



(40) 



(41) 



r exp(^ s /T p0 ) - 1 

-mil ex P i^s/Tpo) 

With t determined by formula (14). It is seen that the temporary behaviour of density 
matrix in Bloch-Redfild approximation do not agreed with above exact solution. But the 
relaxation time of the average energy of the oscillator < n >= J2 n Pn n is coincide with (14) 
This can be shown by multiplying (40) by n and summing over all states: 

1 r 1 

r lexp(u s /T p0 ) - 1 



< n > 



< n > 



(42) 
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FIG. 1. The kinetics of the inverse temperature of the oscillator for various strength of inter- 
action (a,b). The frequency dependance of absorption (c). 
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